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SUMMARY 


A new technique for measuring heat-transfer rates on free-flight 
models in a "ballistic range is described in this report. The accuracy 
of the heat-transfer rates measured in this way is shown to be comparable 
with the accuracy obtained in shock-tube measurements. The specific 
results of the present experiments consist of measurements of the 
stagnation-point heat-transfer rates experienced by a spherical-nosed 
model during flight through air and through carbon dioxide at velocities 
up to 18,000 feet per second. For flight through air these measured 
heat-transfer rates agree well with both the theoretically predicted rates 
and the rates measured in shock tubes. For flight through carbon dioxide 
the heat-transfer rates agree well with the rates measured in a shock tube. 
Two methods of estimating the stagnation-point heat-transfer rates in 
carbon dioxide are compared with the experimental, measurements. At each 
velocity the measured stagnation-point heat-transfer rate in carbon 
dioxide is about the same as the measured heat-transfer rate in air. 


INTRODUCTION 


The purpose of this report is twofold. A method of measuring 
heat-transfer rates in a ballistic range is described, and results are 
presented which were obtained by this method for the stagnation-point 
heat-transfer rate encountered by a spherical-nosed body flying through 
air and through carbon dioxide. The choice of carbon dioxide has been 
prompted by its predicted presence in the atmosphere of Venus. 

Most of the experimental results on heat-transfer rates through 
highly cooled boundary layers have been obtained from measurements made 
in shock tubes. The shock tube correctly simulates the stagnation enthalpy 
which is the parameter of major interest in heat-transfer measurements. 
Unfortunately, the shock tube, as it is ordinarily used, does not correctly 


mulale either the Mach number or the Reynolds number. The ballistic 
nge technique correctly simulates both the Mach number and the stagna- 
on enthalpy. Of course, the actual size of the model is not correctly 
produced in either the range or the shock tube; however, in the range 
is possible to vary the Mach number and Reynolds number independently 
that nonequilibrium effects due to the small size of the model can 
studied . 

All the measurements of stagnation-point heating rates presented in 
is report have been made at an ambient temperature of ^>h0° R and an 
bient pressure of l/lO of an atmosphere. For these conditions the stag- 
tion density is sufficiently high, approximately atmospheric, that the 
ow in the shock layer is in thermal and chemical equilibrium in spite 
the small size of the model. 

The authors would like to acknowledge the assistance of 
. ¥♦ J. Kerwin in the design of the thermocouple model and the coils 
id to pick up the signals transmitted by the model. Interpretation of 
2 signal picked up by the coils required an analysis of the variation 
coil response with the dispersion in flight path from the coil's exact 
iter. The authors wish to acknowledge the contribution of 
W. A. Mersman in carrying out this analysis and the resulting 
put at ions . 


SYMBOLS 


radius of pickup coil 

radius of model coil 

ballistic coefficient 

length of pickup coil 

specific heat at constant pressure 

drag coefficient 

distance between coil centers 

internal energy 

emf induced in pickup coil 

free-energy change 

constant in least squares straight line, defined in equation (A) 
slope of least squares straight line, defined in equation (5) 
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enthalpy 

current In model coil 

equilibrium, constant for constant pressure 
thermal conductivity 
Lewis number 

thickness of copper nose cap 
mutual inductance 
molecular weight 
mass of the model 
number of turns in coil 
Nusselt number 
index of summation 
pressure 

stagnation-point heat-transfer rate 

universal, gas constant 

Reynolds number 

cross-sectional area of model 

temperature 

time 

flow velocity at outer edge of boundary layer behind shock wave, Py 

flight velocity of model 

mole fraction 

distance into copper cap 

distance along copper cap measured from stagnation point 
compre s s ib il ity 

distance along model trajectory 
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a diffusivity of copper 

P velocity gradient measured along the copper cap at stagnation 

point 

(i viscosity 

v kinematic viscosity 

p density 

T dummy variable 

Subscripts 

A pickup coil 

a model coil 

Cu copper cap 

D dissociation 

e external flow 

i component of gas mixture 

0 conditions at z = 0 and t = 0 

P constant pressure 

s stagnation point 

std standard conditions 

t differentiation with respect to time 

w fluid properties at wall 

x differentiation with respect to distance into copper cap 

1 chemical reaction, dissociation of carbon dioxide 

II chemical reaction, dissociation of oxygen 

i conditions ahead of shock wave 


2 


conditions behind shock wave 
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EXPERIMENTAL PROCEDURE 
Description of Apparatus 


All the heat-transfer rates reported herein have been made on small 
models in free flight. The models were launched from a light-gas gun of 
the type described in reference 1. The ballistic range is instrumented 
with spark shadowgraphs and electronic chronographs, as well as the pick-up 
coils described in a subsequent paragraph. 

The construction of the model is shown in figure 1. The model 
consisted of a hemispherical nose with a short cylindrical section and 
a 120° conical tail. The body of the model was made of 2024—T4 aluminum, 
chosen for its high strength-to-weight ratio. The model was held in a 
four-piece nylon sabot during launch. 

A 0.0075-inch-thick copper cap on the nose of the model acted as a 
calorimeter heat gage. A constantan wire, silver soldered to the cap at 
the stagnation point of the model, formed the hot junction of a copper- 
constantan thermocouple. This constantan wire, kept small (0.006-inch 
diameter) near the nose of the model to minimize the conduction of heat 
away from the stagnation point y enlarged to 0.025-inch diameter and 
extended back to the center of the model to form the cold junction of 
the thermocouple with a copper disc. A four-turn coil wound on a boron 
nitride form was connected between the two junctions of the thermocouple 
at the disc in the center of the model and at a point on the periphery 
of the nose cap. Cement was used to join all mating surfaces. To permit 
computation of temperature rise, the electrical resistance of the conplete 
circuit was accurately measured before final assembly. 

The pickup coil was made of B and S gage No. 40 Formvar coated copper 
magnet wire, wound on a Fiberglas form. The coil was a 1-layer, 400 turn, 
center-tapped coil, 4 inches in diameter and l-l/2 inches long. Electro- 
static shielding was incorporated in the coil form to eliminate spurious 
signals resulting from the static charge of the model, the wake, and the 
propellant gases. The inside shielding consisted of B and S gage No. 24 
soft copper wires imbedded in the coil form, spaced l/8 inch apart and 
parallel to the axis of the coil. The shielding wires were 1/32 inch 
from the inside diameter of the coil form. Every other wire was cut off 
at one end of the coil form and connected at the opposite end to a B and 
S gage No. 12 copper wire, which was in the shape of an open-ended copper 
ring, so that there were no shorted turns in the shielding. The elimina- 
tion of shorted turns prevented the model magnetic field from inducing a 
current in the shield which might distort the emf induced in the pickup 
coil. The shielding was completed outside the coil by B and S gage No. 24 
soft copper wires, spaced l/4 inch apart and parallel to the axis of the 
coil. Thus, the entire coil was in an electrostatic Faraday cage. The 
construction of the coil is shown in figure 2. 
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As the stagnation point is heated during the flight of the model, 
the resulting potential difference between the hot junction and the cold 
junction causes a current to flow in the model coil, thus producing a 
magnetic field around the model. This magnetic field in turn induces an 
emf on the pickup coils located along the trajectory of the model. There 
are six pickup coil stations, the first of which is approximately nine 
feet from the muzzle and the last approximately 34 feet from the muzzle. 

Each coil station, located about 3 feet down range of a shadowgraph 
station, consists of a pickup coil connected to a Tektronix type 535 
oscilloscope with a type D high-gain differential input plug-in preampli- 
fier. As the model flies down the range, it interrupts a photobeam which 
triggers the spark gap in the shadowgraph station. When the spark gap 
fires, it simultaneously gives a shadowgraph of the model and triggers 
the sweep of the recording oscilloscope which is connected to the pickup 
coil. The experimental setup is shown schematically in figure 3» The 
magnitude of the voltage induced on a pickup coil is a function of the 
velocity and position of the model with respect to the pickup coil, as 
well as the temperature difference of the thermocouple. Therefore, if 
the velocity and position of the model are known, the temperature differ- 
ence of the thermocouple in the model can be computed from the oscillo- 
scope record. Even though the model coil is contained in a metal body 
and surrounded in flight by an ionized layer of gas, no difficulty in 
transmitting the signal to the pickup coil is experienced. This is because 
the current flowing in the model coil varies slowly over the length of the 
range, so that the frequency of the signal is effectively direct current, 
and the counter electromotiveforces from the currents induced in the body 
shell and the ionized layer are negligible. 

Appendix A contains the details of the computation of the temperature 
history of the copper cap from the oscilloscope records in terms of the 
electrical parameters of the system. It should be mentioned at this point 
that the effective temperature of the hot junction of the thermocouple will 
be the lowest temperature present at the junction, since any currents 
induced by temperature differences at the hot junction are short circuited 
and, hence, produce no net flow through the circuit. In the present case, 
since the constantan stem extended through the copper cap, this lowest 
temperature will be the temperature of the inside or rear surface of 
the copper cap. A typical oscilloscope record is shown in figure 4. The 
disturbances on either side of the signal are caused by the spark discharges 
of the shadowgraph stations. 


Data Reduction 


As indicated in the preceding section, the raw experimental results 
consist of the temperature history of the inside or rear surface of the 
copper cap calorimeter. In order to present the data in the most usable 
form, it is necessary to conpute the heat-transfer rate at the stagnation 
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point from this temperature history. Toward this end the following 
assumptions are made: 

(a) The heat flow is one-dimensional. 

(b) None of the heat transferred to the model penetrates beyond the 
inner surface of the copper cap. 

(c) No temperature gradient exists between the front and the rear 
surfaces of the copper cap. 

The heat-transfer rate is then given by, 

• ^u* dT , 

q = ~d t ( 

where a is the diffusivity of copper, k^ u is the thermal conductivity 
of copper, and l is the thickness of the copper cap. Assumption (a) i 
satisfied since the heat-transfer rate does not vary by more than 5 perc 
over the entire front surface of the copper cap, and assumption (b) is v 
well fulfilled since the boron nitride behind the copper cap is a good h 
insulator. Assumption (c) is intimately related to the finite response 
time of the copper cap calorimeter gage. It is shown in appendix B that 
the error incurred in making assumption (c) is less than 2 percent for 
the present experiments. 

Since there is a certain amount of experimental error in the 
temperatures measured at the various pickup coil stations, the problem 
of fitting the proper temperature history function to the data points 
must be considered. For small changes in the velocity the heat-transfer 
rate may be assumed proportional to the cube of the velocity, 



For flight at sufficiently high Mach numbers the drag coefficient is 
approximately constant, so that 

™ = -BV 2 ('- 

dt u 

where B = (pCpS)/(2mp) is the ballistic coefficient for the model. The 
ballistic coefficient may be determined for each round from a plot of 
l/v versus t, as indicated in reference 2. Equations (l) , (2), and (3) 
may be combined to give, 

T ‘ f + g (0 W 

where 

cc4o 


g = - 


2k Cu 2BV 0 


(5) 
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For each round the temperature is plotted against the square of the 
velocity, and the straight line giving the best least squares fit is found. 

The slope of this line together with the ballistic coefficient, the copper 
cap thickness, and the thermal properties of copper then determine the 
heat- transfer rate q-,. The thermal conductivity and the thermal 
diffusivity of copper were taken from reference 3* 

The procedure outlined above yields a value of q-* corresponding 
to the value of the muzzle velocity, V 0 , for each round fired. Equation 
(2) then gives the value of q which corresponds to a given value of V, 
where V is less than the muzzle velocity, V 0 , but greater than the 
velocity at the last pickup coil. The heating rate q^ is a function of 
the velocity of the model and the temperature of the outer surface of the 
copper cap. Both of these quantities change along the trajectory of the 7 

model. For example, at high muzzle velocities the heating rates will be 
higher, and thus the surface temperature will be higher at each pickup 
coil than it would be for a low muzzle velocity. However, in analyzing 
the data it is desirable to eliminate the surface temperature as a 
variable. Therefore, a specific value of the surface temperature is 
selected, and from each round the values of q and V which correspond 
to the selected surface temperature are determined as follows: The 

difference in temperature between the inner and outer surfaces of the 
copper cap calorimeter is estimated using the graphical technique of 
Schmidt (see ref. L) and the value of % appropriate to the particular 
round under consideration. This temperature difference is then subtracted 
from the specified value of the outer surface wall temperature to give the 
temperature on the inner surface of the copper cap. The velocity which 
corresponds to this temperature is then found from the plot of temperature 
versus the square of the velocity mentioned in the preceding paragraph. 

The proper value of q can then be computed from equation (2) . Thus 
each round fired gives one value of q at the specified wall temperature 
and a particular velocity. For the present series of experiments a wall 
temperature of 900 ° R was selected since this temperature occurred before 
the last pickup coil on the low velocity rounds and after the first 
pickup coil on the high velocity rounds. 


DISCUSSION OF RESULTS 
Heating in Air 


The main purpose of the present experiments was to measure 
stagnation-point heat-transfer rates in carbon dioxide. However, in order 
to compare the results obtained by the present measuring technique with 
the extensive shock tube measurements, several rounds were fired in air 
at a pressure of 76 -mm Hg and a temperature of 5^0° R. The measured 
heat-transfer rates are plotted versus velocity in figure 5* The rates 
predicted by the theory of reference 5 are shown as the solid line , and 
those predicted by the theory of reference 6 are shown as the dashed line. 
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In addition, some of the results of the shock tube measurements reported 
in reference 7> corrected to the present test conditions by the method 
outlined in reference 10, are shown in figure 5* The agreement of the 
present measurements with the theory of reference 5 and the shock- tube 
measurements is good. The rates given "by the theory of reference 6 are 
lower than the measured values. 


The heat-transfer rates are replotted in dimensionless form in fig- 
ure 6 . The heat-transfer parameter is the Nusselt number divided by the 
square root of the Reynolds number. These dimensionless parameters are 
defined by 


Nu 


qy g Pw 

^Sr(^s “ k-w) 


(6a) 


Re 


?£ m Pyf 

\ " v w 


The heat-transfer parameter is then 


Nu 

Tie 


l c p 


V 


^W(hg ky.) 



(6b) 


(7) 


The properties of air necessary for the evaluation of these parameters 
were taken from references 8 and 9* The viscosity was obtained from 
Sutherland’s formula. The velocity gradient at the stagnation point, 0, 
was computed from the Newtonian approximation given in reference 5* It 
is apparent from figure 6 that the present experimental results agree 
very well with the theory of reference 5* Again it is clear that the 
heat-transfer rates predicted by the theory of reference 6 are low. 


Heating in Carbon Dioxide 


Most of the heat-transfer measurements were made in carbon dioxide 
at a pressure of 76-mm Hg and a temperature of 5^0° R. There was a 
maximum of 5 percent air in the test gas because of unavoidable leaks in 
the flight-test chamber* Although the effect of the air on heat-transfer 
rate is unknown, it is very likely quite small for the conditions of 
the present tests. 

The measured heat-transfer rates are plotted versus velocity in 
figure 7. In addition, the experimental results of reference 10, which 
were obtained in a shock tube using a mixture of 90 percent carbon dioxide 
and 10 percent nitrogen, are shown for comparison. These shock tube 
results have been corrected to the ambient conditions of the range 
measurements . 
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No theoretical analysis such as that of reference 5 has been made 
for the case of carbon dioxide* The analysis of reference 5 implicitly 
contains the real gas properties of air in the variation of pp through 
the boundary layer. Since the heat transfer is not sensitive to the 
details of this variation, it is proposed here to use the equations for 
the heat-transfer rate obtained in reference 5 for air, but to substitute 
the thermodynamic and transport properties of carbon dioxide for those of 
air wherever they occur in the final formulas . Although these formulas 
call for the use of properties at the stagnation point, the properties 
just behind the shock wave are used since the difference in the final 
result is negligible. The thermodynamic properties of carbon dioxide at 
the temperatures and pressures encountered behind a normal shock wave at 
velocities up to 20,000 feet per second have been computed as outlined 
in appendix G* The viscosity of carbon dioxide was evaluated from 
Sutherland's formula with C » ^32° R. The properties of carbon dioxide 
at the wall temperature were obtained from reference 9* 

A second method of predicting the theoretical heat-transfer rate 
is given in reference 6. The dashed curve in figure 7 shows the heat- 
transfer rate predicted by this theory. In order to compute the heat- 
transfer rate according to this theory, it is necessary to obtain the 
integral with respect to temperature of the thermal conductivity of 
carbon dioxide from the wall temperature to the stagnation temperatu re. 

These data were obtained from the results of reference 11. 

When compared with the theory of reference the present measure- 
ments give a standard deviation of ±11 percent. The standard deviation 
for the shock-tube measurements for comparison with the same theory is ±12 
percent. Hence, it may be inferred that the accuracy of measurement is 
substantially the same for both systems. 

The heat-transfer rates in carbon dioxide are replotted in 
dimensionless form in figure 8. The dimensionless parameters have been 
determined from equations ( 6 ) and ( 7 )* The stagnation-point velocity 
gradient has been computed from the Newtonian approximation of 
reference 5* 

The theory of reference 6 gives a better estimate than the theory of 
reference 5 for the heating rate in carbon dioxide for Mach numbers from 
12 to 20. However, neither theory satisfactorily predicts the variation 
of heating rate with Mach number. Both theories indicate a decrease in 
heat-transfer parameter with an increase in Mach number, whereas the 
experimental results indicate practically no variation of heat-transfer 
parameter with Mach number; consequently, caution should be exercised in 
using either theory at Mach numbers above 20. 

The effect on the heat-transfer rate of changing the gas from air 
to carbon dioxide may be seen by comparing figures 5 and 7 or figures 
6 and 8. The ambient pressures and temperatures are the same in both cases. 
The inportant result is that the heat-transfer rates measured in carbon 
dioxide are of the same magnitude as those measured in air. 
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CONCLUSIONS 


A novel method of measuring heat-transfer rates in a ballistic 
range has been developed. Trie accuracies attainable with this method 
compare favorably with the accuracies attainable in shock-tube measure- 
ments of heat-transfer rates. 

The potentiality of the system described herein lies in two 
directions. Mach number and Reynolds number can be varied independently 
A and heat-transfer measurements can be made near the base of the model 

^ without fear of interference from a supporting structure. 

2 

T The measured values of stagnation-point heat-transfer rates in air 

agree well with the theory of reference but are higher than the rates 
predicted by the theory of reference 6. The measured values of stagnation- 
point heat-transfer rates in carbon dioxide agree fairly well with the 
theory of reference 6, but lie below the values predicted by the theory 
of reference 5 * Neither theory, however, correctly predicts the nearly 
* constant value of heat-transfer parameter, Nu/VRe , with change in Mach 

number which is shown by the experimental measurements. The measured 
values of heat-transfer rate in carbon dioxide are approximately the 
same as those measured in air. 

Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Jan. 11 , 1961 
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APPENDIX A 

TEMPERATURE HISTORY OP COPPER CAP CALORIMETER GAGE 
EROM COIL RESPONSE 


When the model is in flight and the stagnation point is heated, the 
temperature rises at the hot junction of the thermocouple in the copper 
calorimeter cap- The cold junction at the other end of the constanta n 
rod inside the model far from sources of heat remains effectively at the 
ambient temperature of the model prior to launch during the very short 
time of flight through the range . This temperature difference of the hot 
and cold junctions of the thermocouple results in an emf, causing a 
current to flow in the model coil* The current flow produces a magnetic 
field which induces an emf in the pickup coil as the model passes 
through it. 

The voltage induced in the pickup coil by the moving model is given 
by Faraday f s law 

d(Ml) 
e = - 

dt 

where 

e emf induced in pickup coil 
M mutual inductance between the two coils 


(Al) 


current flowing in the model coil 


If it is assumed that the current in the model coil is constant 
during the short time that it takes the model to pass through the pickup 
coil, the voltage induced in the pickup coil is 

dM 


e = -I 


dt 


(A2) 


From reference 12, the mutual inductance between the two coils may 
be expressed as: 


M = 0.0250 


a 2 N a N A 


D + b 


D - b 


\/(D + b) 2 + A 2 \/ (D - b) 2 + A 2 J 


(A3) 


where 


a radius of smaller coil. In 
A radius of larger coil, in* 
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2b length, of the larger coil, in. 

N a total number of turns on the smaller coil 


total number of turns on the larger coil 
D axial distance between the centers of the coils, in. 


M mutual inductance in microhenries 

Substituting the appropriate values and expressing the mutual 
inductance in terms of the distance between coils yields the equation 


M = 0.3M 


\ 


D +0-75 


D - 0.75 


n/(D + 0.75) 2 + h n/(D - 0.75) 2 + 4 


e = -1.368xlO _S IV 


[(D + 0.75) 2 + ^] 3/2 [(D - 0.75) 2 + ^] 3/2 


(A4) 


Since dM/dt = V(dM/dD), differentiation of equation ( a 4) with 
respect to time followed by substitution into equation (A2) gives 


(A5) 


Because it is more convenient to reduce the data by using the peak 
voltage recorded on the oscilloscope, the above equation is differentiated 
to find the value of D for which e is a maximum. When this value 
of D is substituted in (A5) , 


e max = -0.096X10“ 6 IV (A 6) 

where V is the velocity of the model in inches per second when the 
model is halfway through the pickup coil, I is in amperes, and e re 
is in volts. Equation (A 6) is valid for those cases in which the 
frequency response of the pickup coil and oscilloscope is flat. The flat 
portion of the frequency response can, of course, be extended by properly 
damping the combination of pickup coil inductance and oscilloscope 
capacitance. Although this was not done for the present tests, it would 
be necessary at higher velocities. Preliminary calculations indicate that 
errors from this source will be less than 10 percent for the present 
measurements . 

The peak induced voltage may be measured from the oscilloscope 
record. The model velocity may be computed from the shadowgraph records 
and the chronograph readings. Then the current flowing in the model coil 
may be obtained from equation (A 6) , The product of the calculated current 
and the measured model circuit resistance gives the emf output of the 
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thermocouple . The temperature rise of the copper cap calorimeter may he 
obtained from the thermocouple calibration curve (fig. 9) after the 
necessary correction is made for the cold junction at room temperature. 

The calculations discussed above are valid for on-center shots only. 
Since the launching of a model through the exact center of a coil is but 
a fortuitous occurrence, corrections must be made for off-center shots. 
For off-center shots the voltage induced in the pickup coil is still 
given by equation (A2) , however the mutual inductance between the two 
coils is no longer given by equation (A3)» A more complicated expression 
which involves the lateral distance between the axes of the two coils is 
obtained. This expression has been evaluated on an IBM 650 digital 
computer. 

The correction factor was then experimentallv verified by mutual 
induction measurements as shown by the plotted points in figure 10. The 
correction factor is applied as a reduction in the measured signal 
strength to give the effective signal strength which would have been 
recorded if the model had passed through the coil's center; that is 

eon-center = (correction factor) e 0 ff_ C enter (A7) 

when eoff-center is the measured value of e and eon-center is the 
value to be used in equation (A 6) . 

The magnitude of the signal is also affected by the angle of attack 
of the model. However, the angle of attack usually is slight or even 
undiscemible and is considered to be zero for al 1 the data reported 
herein. 
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APPENDIX B 

TRANSIENT RESPONSE OF CALORIMETER GAGE 


The relationship between the temperature history on the rear face of 
the copper cap and the heat flux into the front face may he determined 
by solving the following boundary-value problem: 


T t (x,t) = dT^Xjt) 
T x (Z,t) = 0 


T(x,0) = 0 


- ^Cu^x(0^) 


do 

(BV 0 t + l) 3 








(Bl) 


where T(x,t) is the temperature at any point x and at any time t , a 
is the diffusivity, kQ u is the coefficient of heat conduction, is 
the heat flux at time t = 0, V Q is the velocity at time t = 0, and B 
is the ballistic coefficient. 


A relationship between T-t(Z,t) and may be found by solving the 
above problem with the Laplace transform method (ref. 13 ) • 

n 2 7t 2 aT 

+ 2 ) (-l) n ' 


T t (Z,t) = 


a 


Lo 


k CuU (BV Q t + l) 3 


n=i 


-3BV 0 do 


0 [ BV 0 ( t - t) + 1]' 


n g 7t g at 

* " l 2 

+ Qo e 

The approximate formula used in the data reduction is 


(B2) 


T t (Z,t) 


a dp 

k CuZ (BV 0 t + l) 3 


(B3) 


The error introduced by the transient effect is thus seen to be given by 
all terms after the first in equation (B2) . 

The relative error is defined as the difference between the exact and 
the approximate heat-transfer rates divided by the exact heat-transfer rate . 
This quantity has been evaluated for V Q = 18,000 feet per second and 



B = 0.003 second -1 which represent the most severe case encountered in 
this series of tests. The results are plotted in figure 11. The error 
from this source over that portion of the range where the measurements 
are made is seen to be less 1.5 percent. 
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APPENDIX C 


IT 


PROPERTIES OF DISSOCIATED CARBON DIOXIDE BEHIND A SHOCK WAVE 


In order to calculate theoretical heat-transfer rates, it is 
necessary to know the properties of the gas "behind the shock wave. Study 
of the magnitudes of the equilibrium constants of possible chemical reac- 
tions at the high temperatures and pressures in the shock layer shows 
that, up to about 6000° K, the major components of a gas consisting ini- 
tially of C0 2 will be C0 2 , CO, 0 2 , and 0. Thus, if the gas is treated 
as a mixture of ideal gases in equilibrium, the required thermodynamic 
properties can be calculated for various temperatures and pressures by 
the same procedure used to calculate the properties of high-temperature 
air (refs. lU and 15)* These properties can then be related to flight 
velocity and ambient density through the normal shock conservation 
equations . 

The first step in this procedure is to determine the proportions of 
the component gases in the mixture . This is done by solving a set of four 
simultaneous equations in four unknowns, consisting of: 

(a) One equation expressing Dalton's law of partial pressures 

(b) One equation of material balance 

(c) Two independent reaction equations 

Dalton's law of partial pressures states that the sum of the partial 
pressures of the components is equal to the total pressure of the mixture. 
In this case, then, it follows that 


*C0 a * %> + *02 + *0 " 1 < cl > 

The material balance equation follows from the fact that in C0 2 
and in any dissociation products derived from it, there will always be 
twice as many oxygen atoms as there are carbon atoms . 


2X C0p + %) + 2X 0 2 + *0 


X C0 2 + ^0 


The two reaction equations are 


co 2 I CO + | 0 2 (C3) 

0 2 t 20 (c4) 
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From the law of mass action it follows that 

(x C0 p)(x 0 p) 1/2 


w 


= K P, 


(C5) 


(Ml! = Kp 

Xq 2 P hi 


(c6) 


where K_ and K are the equilibrium constants for the indicated 
T I T II 

reactions . These constants are calculated for temperatures in the range 
of interest from the ideal gas free energy data of reference 9 according 
to the formula 


*p ■ ex 4 i) (C7) 

When these values of Kp are inserted into equations (C5) and (C6) , 
the equations along with equations (Cl) and (C2) become the set which 
must be solved simultaneously for the mole fractions of the various com- 
ponents as a function of temperature and pressure. By substitution, 
these equations can be reduced to one cubic equation in one variable; if 
this variable is chosen to be, say, the mole fraction of O 2 , then the 
problem is reduced to finding the real root between 0 and 1 of 


2PXq, 


3/2 


V 0 * 172 ' K * 


ft 


0 

(C8) 


The other component mole fractions can then he found by substituting 
this value back into the original equations . 

The next step in the procedure is to use these mole fractions to 
calculate the thermodynamic properties of the mixture as a function of 
temperature and pressure. The properties which are needed to solve the 
normal shock equations are the compressibility, the internal energy, and 
the density. Thus 


Z 2 - 



E 2 


^ Z 2 X ± E ± 
i 


P 2 


Y u 

L- j Z2R/ mQQ 2 T2 


(C9) 


(CIO) 


(cil) 


U W r > 
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where is the internal energy of the ith component of the ideal gas 

obtained from reference 9* 

The final step in the procedure is to relate the properties on 
either side of the shock wave by solving the normal shock conservation 
equations . 

The requirements of conservation of mass, momentum, and energy 
are, respectively, 

PiVi = P 2 V 2 (C12) 

Pi + PiVi 2 = P 2 + P 2 V 2 2 ( 013 ) 

El + + | n 2 - + If + | Vs 2 (014) 

These equations can be combined and simplified according to the 
method outlined in reference 1 6 . Then, when the properties of the 
ambient room-temperature C0 2 ahead of the shock are obtained from 
reference the pressure, density, and temperature behind the shock 
wave can be determined for corresponding values of flight Mach n umb er 
and ambient density as shown in figures 12(a) , 12(b) , and 12(c) . 
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Figure 1 . - Model . 
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Figure 2.- Pickup coil. 
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Figure 3 •- Experimental setup. 
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Amplitude 1 mv/div 
Sweep rate lOOpsec/div 
Launch velocity 11,000 ft/sec 
l/lO atm. CO 2 


Figure 4 .- Typical oscilloscope record. 
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Mach number 

(a) Pressure ratio. 


Figure 12.- Ratios of pressure, density, and temperature across a normal 

shock wave in carbon dioxide. 
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